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Abstract 



The effect of density fluctuations upon light propagation is calculated 
perturbatively in a matter dominated irrotational universe. The starting 
point is the perturbed metric (second order in the perturbation strength), 
while the output is the Hubble diagram. Density fluctuations cause this 
diagram to broaden to a strip. Moreover, the shift of the diagram mimics 
accelerated expansion. 

1 Conventions 

We use the Landau conventions, i.e., we assume H signature for the 

metric. The zeroth component of a four-vector is timelike, the first, second 
and third components are spacelike. Four-vectors are indexed by Latin letters, 
three-vectors by Greek letters. Throughout we use c = 1 units. 

2 The Sachs optical equations 

A light ray may be parametrized along its path by a parameter A so that its 
four-velocity u % is defined by 



It is a null vector, 



u l Ui = 



(2) 



and satisfies the geodesic equation 



u* k u k = . 



(3) 



Let us introduce another independent null vector w l along the path. It satisfies 

vfWi = , (4) 

w 1 ^ = 1 (5) 

and 

w\k^ = . (6) 

Furthermore, let us define two spacelike unit vectors L % a (a, (3 = 1,2) that are 
orthogonal to both the null vectors u l and w % : 

Llm = (7) 

L l a w t = (8) 

L l a Lip = —S a/ 3 (9) 

We also require that the vectors L z a are parallel translated along the light ray: 

Ll, k u k = (10) 

Eqs. @, ([TO} ensure that if Eqs. Q, 0, @, @ are satisfied at one 
single point of the light ray, they will be satisfied at any other point of it as 
well. For the separation of two light rays we assume 

Cm = o (n) 

and 

Cwi = , (12) 

hence we have 

a=l,2 

The coefficients d Q describe the proper separation of the two nearby light rays. 
For their derivative with respect to the path parameter A we have in the most 
general case (since no rotation is possible) 

^ = ^(6^+^)^ (14) 
where w a p is symmetric and traceless: 

= - P ) ■ ^ 
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Now, the geodesic deviation equation 

^ = R^u^e (16) 
leads to the Sachs optical equations for the expansion rate and shear w a p: 

%+® 2 + \ w2 = -\ R i ku3uk 



dw D 



where 



2Qw a p - a, fci L>^ fc 4 (18) 



,2 



W 



(19) 



and Cijki stands for the Weyl tensor. 

The physical situation (e.g., that the two nearby light rays are emitted from 
a point source) is specified by the initial conditions. After having solved these 
equations, governs the change of the cross section A of a light beam: 

f x =2SA. ,20, 

3 Perturbation expansion 

We are going to consider a narrow light beam that propagates in a spacetime 
with a slightly perturbed FRW metric. Hence, the expansion is done over the 
corrections of the metric, the zeroth order being the (flat) FRW metric. We 
expand all the above quantities as follows (if the zeroth order is automatically 
zero, we leave that out): 

W^+^+flg 5 (21) 
Explicitly, according to Ref. [T] 

-a 2 (t)S a/3 

-a 2 (t) (.F5 a p + ^ a (t)t 2 Q F a ,p) (22) 











= O.flS 




(2) 

9oo 


= 0,g% 


= 


(2) 
9a/3 


= -« 2 w 





81 
2800 



a 2 (t)t 4 Q [l9F;JF 7!(3 - 12F^F n F a ^ + 3 ((F<;) 2 - F;JF;^ S a0 



Here to stands for the present time, while a(t) = (t/to) 2 ^ 3 denotes the standard 
Friedmann-Robertson- Walker scale factor in a matter dominated homogeneous, 
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isotropic, flat universe. Indices are raised and lowered by the Eucledian metric 
5 a p. Our consideration is restricted to scalar perturbations and the expanding 
mode. It is characterized by the scalar function F({x a }) which depends only 
on the spatial coordinates. 

r^^+T^+Tf (23) 

R ij =B§ ) +B$ ) +B$ ) (24) 

Ciju = CjjM + C\) kl (25) 

u l = u 1 ^ + u 1 ^ + u l( V (26) 

w i = w i(0) +w iW +w iW (27) 

^ = L%°) + LjW + V^> (28) 
6 = (O) + Q(i) + @(2) (29) 

WaP=wW+wW (30) 

A = A^+A^+A^ (31) 

The trajectory of the light ray will be given parametrically as (t(X), r(X),'&(X),<p(X)), 
which will also be expanded as 



t(X) = 


--t^(X)+t^(X)+t^(X) . 


(32) 


r(A) = 


r (o) (A)+r (D (A) + r (2) (A) _ 


(33) 




= i?W+i?W(A) + ^ 2 )(A) . 


(34) 


4>(X) 


= 0(°)+^ 1 )(A) + 0( 2 )(A) . 


(35) 



Note that we shall use throughout a comoving synchronous gauge, since we 
consider a dust-filled universe. Hence, 500 = 1 and goa — to all orders (using 
c=l). 

As we consider the perturbation of a light ray, all quantities depend on the 
parameter A (we shall use mainly f(°) instead). In case of the field variables this 
means explicitly that e.g. 

9ij = 9ij{x J (X)) , (36) 

i.e., we take the arguments along the light ray. This dependence should also be 
taken account when performing the perturbation expansion. 
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4 Zeroth order 

In zeroth order we write for the line element 

ds 2 = dt 2 - a 2 (t) (dr 2 + r 2 (dd 2 + sin 2 M<i> 2 ) ) (37) 

thus is diagonal, and 

9^ = 1 (38) 

9^=-a\t) (39) 

g<® = -a 2 (t)r 2 (40) 

C = -« 2 (t)r 2 sin 2 (41) 

In a matter dominated universe with zero pressure we have 

2 

aW = a °(^) 3 (42) 

where ao stands for the scale factor at the present age to and r is a parameter 
depending on the previous history of the universe. In what follows, we set r = 0, 
which means that f is not measured from the Big Bang, and thus t is not the 
actual age of the universe. Certainly, one may change in the final expressions t 

to t — t again. Also, we take ao = 1. The nonzero components of the Christoffcl 
symbol are 

lf r °) = ad (43) 

r#> = r 2 ad (44) 

rj° } = r 2 sin 2 tiaa (45) 

r r(0) _ r r(0) _ pi?(0) _ r f (0) _ r 0(O) _ r <«0) _ a afi x 

= -r (47) 

r;f=-r S m 2 ti (48) 

ry = -sintfeostf (49) 

ptf(O) _ r 0(O) _ p0(O) _ r 0(O) _ 1 

rft 0) =ra o) = cottf (5i) 



#0 

The Ricci tensor is diagonal with the elements 



R?t> = -3- (52) 

-ai + 2a 2 (53) 

4°j = r 2 (aa + 2a 2 ) (54) 

= r 2 sin 2 (ad + 2a 2 ) (55) 
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The Weyl tensor vanishes. The results for a light beam emitted from a point 
source at time t\ and arriving at the origin at to are the following: 

u*<°> = l a (56) 



« r(0) =-\ (57) 



1 



a 



u 



m = u m = o (58) 



For the trajectory we get 



r (0) = -l (59) 



or 



r(°>(t) = - f tdt =3(^-4^) (W'l 
J t0 a \a(t ) a(t) ' 



while and are constant. We also introduce the radial coordinate dis- 
tance from the source: 

*) = rm < t ')- m « = 3 (4-«)- < 61 » 

For the expansion rate we have 

e<°) = 4 + 4- = M (62) 

a z a z d a l a 
and thus the cross section area of the beam is 

A™ = na 2 d 2 (63) 

The intensity of the light beam is then 

j(Q)- £ " "° ( ° )2 fa) L a ^ f64) 

47ru ( ) 2 (ti)A(°)(io) 47ra 4 (t )rf 2 (io) 

where L stands for the absolute luminosity (radiated energy pro time unit). It 
is customary to express it in terms of the red shift 

measured wavelength — emitted wavelength u°^(ti) ^ a(t ) wgr\ 
emitted wavelength u°(°)(i ) 

Since the scale factor a is proportional to t 2 ' 3 , we have 

L H 2 1 L 2 / 1 1,7 13 



^ 4 (i + z )(yr + ^-l) 2 4tt 2z 16 32 
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or 



log/(°) = - log (^f\ + 2 logz + \z - ^ + ^ + ... (67) 



The Hubble parameter H is given by 

H — - 
a 



2 



(68) 



Note again that to is not the actual age of the universe. 

We shall also need the zeroth order expression of w l and L % a . They are given 

by 



„*(o) - a 



2 
1 

2 

w m = w m = (71) 



W H0) 



(69) 
(70) 



and 



L*g> = L[f = (72) 

Lf (0) = (74) 

L* (0) = (75) 

r^(0) _ 1 ( 7 a) 

2 ~ arWsin^W 1 > 



5 First order 



As w a p vanishes in zeroth order, Eq.(17) implies that for a second order calcu- 
lation of O one needs the first order corrections of w a p and Cijki- In contrast, 
corrections of Rjk and u 3 are needed up to second order. As the corrections to 
the metric are given, we consider it and all of its derivatives as known quantities 
and seek for corrections of the path, expansion and shear of the light beam. In 
this section we consider the first order corrections. As the Sachs optical equa- 
tions refer to a certain light trajectory, the deformation of that trajectory should 
also be taken into account. Therefore, we start with the equations Q-pl). For 



the first order correction of u 3 we have 



d9^ 



2g^^uf } + fl £V«VC°) + -^r a^C W) = . (77) 

J J J (jx 

Inserting the zeroth order expressions we have 

2dt^ dr« g$ 2d m „ 

-^+ 2 ^7^ + ^-V () =0. 78 
a dX dX a 4 a 3 
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Similarly, the geodesic equation implies 



du^ 1 



dX JK Jfc 



ox 1 



Inserting the zeroth order expressions this leads to 



d 2 t^ hdrW 



dX 



dX 2 



dX 2 
a dtW 



dA 1 



(V 



dX 



- 2 



dX 
dX 2 
dX 2 



(1) _ 



,(ar'(°)) fflW 
J a 2 r(°) dX 

' a 2 r(°) dA 



-rJ<V<°>u*<°> 



(79) 

(80) 
(81) 
(82) 
(83) 



As the right hand sides of these equation s a re in prin ci ple known functions 
of A (or, equivalently, of C = i (0) , cf. Eq.g) Eqs.g, @ can be directly 
integrated to get 



Mi) 



(0 



1 



Q 2 r (0)2 



dC a V°> a rJ<V< >« fc <°> 



(84) 



Here the constant of integration (or, equivalently, the lower limit of the integral) 
is determined by the criterion that ^^'(Q must not diverge at £ = t . A 
divergence corresponds to a zeroth order change in the direction, describing a 
light ray not hitting the observation point (i.e., the origin). As Eq.(84) shows, 
such divergence does not occur, and the direction of the light ray is not even 
changed at £ = to- It does change (due to the perturbation), however, at the 
source, i.e., at £ = t\. 

A further integration yields the first correction to the angle i?, namely, 
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_J f 1 a 3 r (0)2 r *(l).J(0) k(0) 



(85) 



+ 



The integration constant has been determined by the condition that #W = 
as the initial position is fixed. Similarly we have for the other angle 



(1) (C) 



1 



a 2 r (0)2 



(86) 



1 



r(°)(C) 
1 

r( n )(t!) 



to 



• dC > a 3 r (0)2 r *(l) u i(0) u fc(0) 



(87) 



+ 



ti 



All quantities above are considered as functions of C = . 

Eqs. ( 78 1 , ( 80 ) , ( 81 ) are not independent, hence, it is advisable to use Ens. ( 78 1 
and (80). Then, expressing < ^^- from Eq.(78) and putting it into Eq.(80) we 
have 



i?D + Aid) 



(i) _ 



_ r t(1) 7/j(°)7; fc (°) 
5 i/rr 1 jk " " 



Again, a dot means derivative with respect to ( = t^°\ The solution of Eq.(88 1 
can be written as 



t 



(i) 



t 



c 



where 



f(0 = -~9$-a 2 rW u m,M°) 



3 i i9 2 F 
10 oC dr 2 



(89) 



(90) 



while a and /? are integration constants. Note that in Eq.(90l the arguments of 
d 2 F/dr 2 are to be taken along the light ray, hence, it also depends on (. We 
also made use of the fact that a oc £ 2 / 3 . Now we have 

eftW 1 dtW 



a dC, 



13 
a \ 5 



(91) 
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n ' m =-^{l £ d( m + l ci £ * f!/(e) + i + 5° - j <92) 

and 

r (1) =-|^ [ C dU(0 + l C — f>^ f /(0- / C ^^|--C-^C-+7(93) 
5 a 7 to 5 a J t J t 2a a a 

where 7 is a further integration constant. The three integration constants a, ft, 
7 are determined by the conditions 

i (1) (*o) = (94) 
r {1) (to) = (95) 
£ (1) (ti) = (96) 



We get 



°- iAM 1 -(£)') /(a ,97) 

" - -iftr*( i -(s) , ) /t<) <98) 



(99) 



The first order correction of the expansion rate satisfies the equation (cf. 
Eq.©) 

dA 2 Jfc J * 2 <9x' 

Inserting the zeroth order terms on the left hand side, we have 

i^ + 2 MeC) = x «(c, (101, 

a aQ a z a 



where X^KO stands for the right hand side of Eq. ( 100 ) , that is a known function 
by now. Explicitly, it is 

_ s 

X W (C) = 3 AF ( 102 ) 



10 V*o 

2i_/^\- 3 ^_4i_/^\-^ F 4i_/^y 5 r< n 

2 1 fc\- 5 F( tl )-F(t ) 41 ^c^~ 5 fe) 3 r^fi'j r 
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The solution of Eq.(lOl) can be readily written as 

0(1) = ^/> a3d2 X (1) (O (103) 

Note that the integration constant (i.e., the lower integration border) is fixed by 
the condition that no singularity can appear at £ = t\ since that would amount 
to a change of the light source position. The first order correction of the beam 
area satisfies the equation (cf. Eq.(20)) 

^ = 2e(°> J 4« + 26( 1 ) J 4(°' (104) 
or, inserting the explicit zeroth order expressions 

I^_ 2 M A (i) =2e (^a 2 d 2 . (105) 
a dQ a z d 

Its solution is 



A {1) = 2Qa 2 d 2 £ d£ a9 (1) = 2Qa 2 d 2 ^-^ £ d£ 



a 3 d 2 X (1 H0 



C d^a 3 d X {1) (0j ■ (106) 

The integration constant is determined by the condition that the solid angle at 
small distances is given, i.e., A^/(a 2 d 2 ) vanishes when £ — > t\. Explicitly we 
have 

AW = A^ltS — _J V _5 J -(F(h)-F{to)) 



t 3 — t 



t 3 t 



It-^-t / deAF(o(-i+(4+*f)r*-4*?r s ) [ • (107) 



Let us sketch how to get the corrections to the Hubble diagram. We have now 
expressions for A and it* parametrized by t\. (Note that to is fixed.) The 
intensity and redshift are expressed as 



■47T („t( tl ))2 A(f|l) 

and 



(109) 
«'(*o) 
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respectively. Hence, the Hubble diagram I(z) is given in a parametric form. 

In the second order calculation we shall also need the first order correction 
to w a p. From Eq.(18) we get 



dw 



(i) 



(110) 



i.e. 



(i) 



^ + 2M^) =Ca4 (0)^(0) M M0)4(0) 



(111) 



The solution sounds 



1 



w 



(i) _ 



a*d 2 cl]lLf\W u m L f) 



(112) 



Again, the integration constant is fixed by the condition that no singularity may 
appear at £ = t\. 

For completeness, below we present the explicit expressions for all the rele- 
vant first order corrections. 



t«(C) 



,(1) 



(0 



— C ^(0 
10 s vs; 



lo 1 " 1 

3 
5 



to 

< 



di 



to 



to 



di 



"3 +3 



C 



'n 



3 
1 

3 



F(0~ 



io* n,i 



-3 0- 

4 - 



4.3 _f3 



di 



to 



(113) 



f(0 



(114) 



+ 



20 
to 



. 1 n 



dF 
dr 



(to) 



dF 
dr 



(0 
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to F(t ) 



to 



F(0 



3 qF(t ) - 0f(o 3 qF(t ) - qF(h) 



10 



10 



tf«(C) 



3 1 



5 t 



1- ^ 



(4)' 



.(o)( C ) 



Co 



to 



9 t§ 



3 1 



5 t 



1 - a- 



ft) 



to 



Q2p r (0)2 ^3^ \ 

20 \jM~ ry ' drdti + d 2 rdtf J 



„(o) 



3 to a 



5r(°) V*o 



drdfi 



) 



dF 



12 



1 



r(°)(ti 

3 Jo_ 
5r(°) 



/2 
c 



20 r(°) 2 



3 t 

5r (0)2 



3.F 



,(o) 



dF 

' 

d 2 F \ 

dr&d) 



AO) 



d 2 F 

drdd 



r (0)2 Q3p \ 

^TWrlM] 



3 <9F 



ft? 



drd-d 



+ 



,(0)2 ^3^ \ 

~2~~ 5 2 r<9t9 J 



at? 



<9r<9i9 / 



1 



r(°)(C) sin 2 ^y t0 

i 

to 



3 Jo_ 
5r(°) 



l 



20 r(°) 2 

)2 



<9F 

90 



1 1 

2^ 

_-.(0) 



i 
to 

d 2 F 
drdcf) 



3 <9F 

r (0)2 \ 

"T" <9 2 r<90 J 



1 



r(°)(ti) sin 2 



r (0)^n_ 
9 t 2 , /<9F 



<9F 



3 t 
5r(°) 



— / 

sm 2 1? it, 



3 to 



5 r(°) 2 
3 

To to 

1 

10 

2 1 

5to~ 

9 2 
20 to 

1 

To 



to 



dF 

20r(°) 3 
<9F 



(0 



20r(°) 2 \d(f> 

2p\ X 



.(0) 



r (0)2 Q3 F \ 



+ 



_ r (0)^ 



dr<90 ' 2 



d 2 rd(j> J 



2 Vt 



<9F 



,(0)^ + 



90 

r (0)2 Q3 F \ 



dr<90 2 d 2 rd(f) J 



r io)^l\ 
drdcf) J 



1 1 

2^°I 



<9F 



<9F 
<9r 



3 

To 



^(0 



2t s + 3C 



C 



to 

d 2 F 



to 



t 



§ F(<l) + To 



2tf + 3C* 



F(t ) 



F(0 + 



(0 



dr 2 v,y 10 
ti\ f 4t f +C f 



to 



1 1 

5to" 

_1 <9F 
dr 



to 

(0 



2tn 3 



3C S 



to 



J. 3 _ J- 3 

t t l 



^(tl) 



41 

5 t V*o, 



/ ^ 

J t 



i 

To 
i l 

5t^ 



4t|_+C3 

|3 _ J- 3 
L L l 



nto) 



Jtn 



t 







ti 



.*0 



13 



/ (1) (C) = 



+ 



r(°) 2 V^o 
3 t (Z 



i*W(C) = 



5r(°) V*o 



1 



to 

OF 



9 4 (dF 



,(o) 



20 r(°) 2 V dd ' drdti 
d 2 F_\ _ 1 / ; 



8 2 F r(°) 2 d 3 F \ 



2 5 2 r<9i9 J 



,(o) 



+ 



c 

r (O) 2 sin 2 i? V*o 



C 



9 t 2 



dF 



dF 



(119) 



wii(C) 



3 Jo 
5r(°) V*o 

-w 22 (C) 



<9f 

00 



,(o) 



20r(°) 2 
d 2 F\ 1 

to 



_ r (0). 



a 2 F r(°) 2 a 3 F ^ 



drd(j) 

dF 
~d~4> 



2 a 2 r9</> y 



(120) 
(121) 



30 C 2 sin 2 i? 
w\i{Q = W2l(C) 



1 



15 C 2 sin 2 




a 2 ^ . 2 a a 2 F . a bf 

sin 17 -^tt + smtfcosw 



<9i9 2 



00 



(122) 



. d 2 F 
cost)— — sini? 



>-(*)' 

The first order corrections to the Hubble diagram are given by 



dF 



+ 



i 



In/ - ln/ (0) (z) 



10 V*o 



<9r 

(Ffa) - F(to)) 



<9F 



(to) 



(123) 



10 



*o 5 f£(*o)-*?^(*i) 



3 dF i 



+ 




/ _ 1 _ 1 \ / _ 2 _ 2 \ / _ 2 _ 1 _ 1 _ 2 \ > 



f 3 — t 



{F(h) - F(t )) 



*° 2 2 



r 1 1 

D i3 J3 

l 4 



ti 



(124) 
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or, in terms of the luminosity distance, 



AW 
A(o) 



= rto 1 - 



3 

10*° 



_ 2 I 
J. 3i3 



/ 3 — /• 



ti 3 *o 3 ) ( ^i 3 *0 ' 



J. 3 _ J- 3 

_ 2 _ 1 _ 1 _ 2 

2i^ -f~ i-^ tg "I - 2£q 



/ 3 / 
2 I 



/•to 



"a ''O 

9 i 



i 3 



tj 5 J ti 



t 3 — t 3 



<) f to / 11,11 

/ ^AF(e)(-i+(t|+tf)r s -^*ir 

t c 1 1 



j(0) 



10° + 



1 + z 



- 2- 



z (2 + yT+j + 2(1 + z)) 



io^vrr^-1 "i + yrr^+(i + z) + (yrr^) 3 + (i + z) 2 

« (1 + zf /2 rto 



(F(h) - F(to)) 



2f 



(1 + Z) 5/2 - 1 Ao(l+2)- 3 / 2 



9 yrr^ 
5* ^r+^-i 



to(l+ 2 )- 3 / 2 



«AF(0 - 1 + *o (1 



1 

VTTz~ 



(F(ti) - F(t )) 



to 

to(l+z)- 3 / 2 



r*- 



2 



vT 



(125) 



Note that Eg.(|123[) is just the Sachs- Wolfe effect 0. 



6 Second order 

First of all, we need the second order correction to the path and the four velocity 
u\ Again, we use both Eq.Q and Eq.|3]). We get from Eq.Q 

2 drW _ 2a p, = _s| _ 2^(1)^(0) _ A^Vd) 

a dA dA a 3 a 4 3iJ ^ J 

floC ) 1 a (1) 1 Q2 (0) 

dx l a 4 9x z 2a 4 d^dx™ v 1 

Similarly, Eq.([3]) leads to 

dX 2 a dX V a 3 a 4 J Jfc 
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(128) 



n r t(0) 

-2r*<V«u fc < > - r*4 0) « i(1) «* (1) ~ 2-^rx l{1) u^)u k ^ 

J* 3 It 5a;' 

n F t(l) -, o2 r *(0) 

^(1)^(0)^(0) _ i^*_ a .l(l) a .m(l) uJ -(0) u fc(0) (12?) 

&r ; 2dx l dx m v 7 

dA 2 a 3 dA a 2 dA U 4 « 5 / ~ 3 ' fc 

sr r(1) 1 a 2 r r(0) 

^(1)^(0)^(0) _ ilii*_ a .l(l) a .m(l) uJ -(0) u fc(0) 

ox 2 ox l ax m 

<P¥» {aM)d*&> _ ^(2)^(0)^(0) 

o r i?(0) 

-2rfiV«u fe (°) - rffV^VW - 2^-z« 1 V< 1 >u*< > 

Jft Jfc Q x l 

_^^_ x HD^(o) u Mo) _ 1^1^^(1)^(1)^(0)^(0) (129) 

dA 2 +%2 r (0) dA 
a r 0(l) , D2 r 0(O) 

_^_^(i) u i(o) n fe(o) „ 1^^^(1)^(1)^(0)^(0) (130) 
dx l 2dx l dx m K ' 

Since the left hand sides are of the same form as in the case of the first cor- 
rections, while the right hand sides are already known, these equations can be 
solved by the same technique as Eqs.(78l-(83|. Once the path and the four 
velocity are known to second order, one can determine the second correction to 
the expansion rate as well. The Sachs optical equation implies (cf. Eq.(lOl)) 

ldQ™ (ad) m / n \ . , , . 



with 



1 1 rlFf^' 

x (2) = _^ R f) u m u m _ R w u m u m _ l --^x l ^ u ^u k ^ 

jfc 2 ■ Jfc dx l 

_l°p* x im u J(o) u Ho) _ 11^^^(1)^(1)^(0)^(0) 
2 dx l 4dx l dx m 

_l(i) 2 _ e (i) 2 . (132) 
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As before, the solution of Eq.( 131 ) is 

8 (2) = ^jf #<W a) (0- (133) 
For the second order correction of the beam area we get the equation 

IdAV _ 2 M A (2) = 2 Q(2) na 2 d 2 + 2 Q(1) A (X) (134) 

a dQ a z d 
Its solution is 

AW = 2a 2 d 2 £ di (a09< 2 ) + -^6 (1) A (1) ^ . (135) 
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7 How does the Hubble diagram change? 



We consider now the question how the Hubble diagram changes due to the 
perturbation. As the Hubble diagram is given in a parametric form, we calculate 
the difference 



ln/-ln/ (0) (z) , 



(136) 



i.e., we compare the perturbed intensity with that of the nonperturbed one, 
taking the latter at the perturbed redshift value. Namely, 



z = 2 W + zW + 



where 



zW= a (t )[ U t W(t 1 )-(l + z(°)yW(io) , 
z (2) =a(i ) [^(^-(l W°V (2) (io)" 
-a\t )u^\t ) [u*M(h) (1 + z«V (1) (*o)' . 
In this way we obtain 

z (i) (l + zC))^ 1 ) 



(137) 

(138) 
(139) 

(140) 



ln/-ln/ (0) (z) = 



(i + z(o))(vTT^r-i) mlfViT^-ij 



,(2) 



(i + z(o))(vr+^)-i) 9ti(vT+z^-i)' 

(2-3v / l + z(°))z( 1 ) 2 (l + z W)2 A (i)2 



4(l + z(°)) 2 (yi + zM - l) 162^ (vT+zW - l)' 
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A Calculating A® 



(141) 



First we make the expression of x^ 1 ' more explicite. According to Eqs.(lOO) 



(101) 



I , qd(0) 

x«(C) = -\n^u m ^ m - R?M l) u m ~ J-Sr^ (1)uJ ' (0)ufe(0) ( 142 ) 



2 <9x' 
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Inserting the zeroth order expressions, we have 



X 
Here 



2a 4 ' 



a 

t2 ' 



ad + 2d 



r(l) 



d(1) 
•Kit 



— AS — 30— -B 
10 a 



a a — 4da 



(144) 



i (1 ^143) 



fl£> = -2(a a) 1 / 2 



as 

9r 



(145) 



rr 10 a ^ ' 5 a b dr 2 



12< "'" 1 288 (— )V 



-< 2 



B 



Combining Eqs.([l43|-p46|), (p9|-(|93|) we get 

x (1) (C) = 



1 -Wx _L»(i)_ — 



2a 2 " + a 3^ 2a 4 



4 1 



1 s $ 20 



3a 2 C 11/3 7 t0 d ^ 3f ^ 3a 4 C 2 ' 9 a 2 C n / 3 



(146) 



(147) 
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